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Abstract. In this paper we consider the eigenvalue problem for piezoelectric shallow 
shells and we show that, as the thickness of the shell goes to zero, the eigensolutions 
of the three-dimensional piezoelectric shells converge to the eigensolutions of a two- 
dimensional eigenvalue problem. 
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1. Introduction 

Lower dimensional models of shells are preferred in numerical computations to three- 
dimensional models when the thickness of the shells is 'very small'. A lot of work has 
been done on the lower dimensional approximation of boundary value and eigenvalue 
problem for elastic plates and shells (cf. |2 3 4 5 6 8 9|). Recently some work has been 
done on the lower dimensional approximation of boundary value problem for piezoelec- 
tric shells (cf. Q). 

In this paper, we would like to study the limiting behaviour of the eigenvalue problems 
for thin piezoelectric shallow shells. We begin with a brief description of the problem and 
describe the results obtained. 

Let & £ = 4> £ (£2 £ ),£2 £ = co x (-£,£) with CO C K 2 , and the mapping 4> £ : H £ — > R 3 is 
given by 



4> £ (x £ ) = (xi,x 2 ,e9{x u x 2 )) +x%a%(xi,x 2 ) 

for all x £ = (xi,X2,x%) € Q. £ , where 9 is an injective mapping of class C 3 and a £ is a unit 
normal vector to the middle surface <t> £ (to) of the shell. Let 70, y e C dco with meas(yo) > 
and meas(7 e ) > 0. Let = 4> £ (yo x (—£,£)) and let f £ = 4> £ (y e x (—£,£)), The shell is 
clamped along the portion of the lateral surface. 

Then the variational form of the eigenvalue problem consists of finding the displace- 
ment vector m £ , the electric potential <p £ and <f; £ € R satisfying eq. \2.2\\ . We then show 
that the component of the eigenvector involving the electric potential <p £ can be uniquely 
determined in terms of the displacement vector u £ and the problem thus reduces to finding 
(m £ ,^ £ ) satisfying equations J2.43t and ( I2.44t . 

After making appropriate scalings on the data and the unknowns, we transfer the prob- 
lem to a domain £1 = cox (—1,1) which is independent of £. Then we show that the scaled 
eigensolutions converge to the solutions of a two-dimensional eigenvalue problem J6.50l >. 
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2. The three-dimensional problem 

Throughout this paper, Latin indices vary over the set {1,2,3} and Greek indices over 
the set {1,2} for the components of vectors and tensors. The summation over repeated 
indices will be used. 

Let ft) C R 2 be a bounded domain with a Lipschitz continuous boundary y and let CO 
lie locally on one side of y. Let Yo,Ye C dco with meas(yo) > and meas(y.) > 0. Let 
7i = <9©\yo and y s = dco\y e . For each e > 0, we define the sets 

Q. £ = co x (—£,£), r ±,£ = co x {±e}, Tq = 7o x (-e,e), 
r i = Yi x (-£,£), r £ = y,x (-£,£), rf = y s x (-£,£). 

Let x £ = (x\ ,X2,xf) be a generic point on £2 £ and let d a = <3„ = and <9| = 

We assume that for each £, we are given a function 9 £ : co — > K of class C 3 . We then 
define the map £ : co — > M 3 by 

^ £ (xi,x 2 ) = (xi,x 2 ,0 £ (xi,x 2 )) for all (xi,x 2 ) € 0). (2.1) 

At each point of the surface S £ — (j) £ (co), we define the normal vector 

a £ = (\d l e £ \ 2 + \d 2 e £ \ 2 + i)- l / 2 (-d l e £ ,-d 2 e £ ,i). 

For each £ > 0, we define the mapping 4> £ : i2 £ — > R 3 by 

4> £ (jc £ ) = £ (xi,x 2 )+xfa £ (xi,x 2 ) forall^ e eO e . (2.2) 

It can be shown that there exists an £o > such that the mappings <t> £ : £2 £ — > <I> £ (i2 £ ) 
are C 1 diffeomorphisms for all < £ < £o- The set Cl £ — 4> £ (i2 £ ) is the reference config- 
uration of the shell. For < £ < £o, we define the sets 

^ = *«(r t .«), f £ = ^> £ (r £ ), ff = 4>(rf), f| / = ffuf ±e ! 
f e £ = 4>(rf), ff = 4>(r £ ), ff D = ffuf ±£ 

and we define vectors gf and g , £ by the relations 
gf = d £ <t> £ and gJ' E -gf = 8i 

which form the covariant and contravariant basis respectively of the tangent plane of 
4> £ (i2 £ ) at 4> £ (x £ ). The covariant and contravariant metric tensors are given respectively 
by 

gfj = gf-g e j and g^ = g^.gj>e. 
The Christoffel symbols are defined by 

r?f = g p ' e -df g f. 

Note however that when the set Q. £ is of the special form Q. £ = CO x (—£,£) and the 
mapping <t> £ is of the form 12. 2\ . the following relations hold: 
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The volume element is given by y/g £ dx £ where 
S £ = det(s?.). 

It can be shown that there exist constants g\ and g2 such that 

< gi < g e < g 2 



(2.3) 



for < e < £o. 

Let A'j kl ' e ,P' jk ' £ and S l ^ £ be the elastic, piezoelectric and dielectric tensors respec- 
tively. We assume that the material of the shell is homogeneous and isotropic. Then the 
elasticity tensor is given by 



(2.4) 



where X and jj, are the Lame constants of the material. 

These tensors satisfy the following coercive relations. There exists a constant C > 
such that for all symmetric tensors (My) and for any vector (f,) £ R 3 , 



i kl '%ti >cttj- 

.7=1 

Moreover we have the symmetries 

j^ijkl,e _j^klij,e _j^jikl,e gkl,z _ ^kl,e pijkx _ pkij,e 

Then the eigenvalue problem consists of finding (u £ , <p £ , | £ ) such that 
-divd £ (i2 £ ,<p £ ) = E, £ u £ in Cl £ 



a £ (u £ ,(p £ )v = 0onr; 
u £ = on Tq 



div£ £ (w £ ,<p £ )=0in 6 £ 
D £ (i2 £ ,<p £ )v = 0onf £ 
r = 0onf £ D . 



where 



where e £ .(ii £ ) = \(dfu e j + dfuf),df = d/dx £ <mdE k (<p £ ) = -v(<f> £ )- 
We define the spaces 

y £ = {ve(// 1 (n £ )) 3 ,v| f§ =0}, 
❖ £ = {v/e// 1 (i2 £ ),vir| fE =0}. 



(2.5) 
(2.6) 



(2.7) 



(2.8) 



(2.9) 
(2.10) 



(2.11) 
(2.12) 
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Then the variational form of systems j2.7t and J2.8t i is to find (u e ,(f> £ ,E, £ ) £ V £ x *P £ x 
R such that 

a £ ((t5 £ ,^ £ ),(y £ ,V/ £ )) = <^ £ /" £ (y £ ,^) for all (v £ , y/ £ ) <E V £ x 4"\ (2.13) 



where 



'-((u e ,p),(v e ,r))= i A'^' £ e 1 / (« £ )4(v £ )di £ 



(2.14) 

^(v 8 ,^) = / « £ -v £ di £ . (2.15) 

£ — £ . 

Since the mappings <t> £ : £2 — > £2 are assumed to be C diffeomorphisms, the corre- 
spondences that associate with every element v £ £ V £ , the vector 

v £ = v £ -4> £ :£2 £ ^R 3 

and with every element \j/ £ 6 *P e , the function 

y/ £ = ^ • 4> £ : £2 £ -> R 

induce bijections between the spaces V £ and V £ , and the spaces *P £ and l P £ respectively, 
where 

y £ = {v £ e (// 1 (n £ )) 3 |v £ = 0onrg} ; (2.16) 
^ £ = {y/ £ eH 1 ^ 6 )!^ = 0onr £ £) }. (2.17) 

Then we have 

d £ v £ (x £ ) = (d £ v £ )(g i - £ ) J , (2.18) 

i iJ m#)=4\\i(y e )bf*)i(g l *)j, (2.i9) 

where 

^(v £ ) = i(^ + 5/vf)-r^v £ . (2.20) 

Then the variational form (12. 1 3i posed on the domain £2 £ is to find (u £ ,(p £ ,^ £ ) EV £ x 
f E xl such that 



where 



a £ ((u £ ,(p £ ),(v £ ,y/ £ )) = ^ £ l £ {v £ ,\i/ £ ) forall(v £ ,v/ £ )ey £ x l P £ , (2.21) 



; (( M £ ,<p £ ),(v £ ,V/ £ )) = ^A^ / ' £ e £ ||/ (v £ )4 7 (v £ ) v ^dx £ 
+ / <r'^ £ (p £ d £ ^ £ v/^dx £ 
+ j a / mi *{d £ m <p £ e £ ^v £ ) 

-d> £ 4 ; ( M £ ))V?dx £ , (2.22) 
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A«"' e =A ijkl ' e (g p ' s )i ■ (g«> e )j ■ (g r - e ) k ■ (g s - £ ) h 
P Pqm =P ijk ' e (g P ^)i-(g q ' e )j-(g r ' e )k- 



(2.23) 

(2.24) 
(2.25) 
(2.26) 



Using the relations A2.3I . ( 12. 5> and ( 12. 6> . it can be shown that there exists a constant 
C > such that for all symmetric tensor (Mi A and for any vector (f,-) S M 3 , 



A^MuMtj > C £ (M !; ) 2 

3 

g ihe titj>cY,ti- 



(2.27) 



(2.28) 



Clearly the bilinear form associated with the left-hand side of (I2.2H is elliptic. Hence 
by Lax-Milgram theorem, given f £ € V' £ and /i £ G V P' £ , there exists a unique (m £ , <p £ ) 6 
V £ x ^ such that 

a £ (( M £ ,<p £ ),(v £ ,v/ £ )) = ((f,h £ ),{v e ,\^)) VV £ x l P £ ey £ x l I' £ . (2.29) 

In particular, for each / £ e (L 2 (£2 £ )) 3 , there exists a unique solution (m £ , <p £ ) € V £ x 
such that 

a £ (( M £ ,<p £ ),(v £ ,v/ £ )) = / / £ v £ v /^ £ dx £ Vv £ x v/ £ eV £ x>I /£ . (2.30) 
This is equivalent to the following equations. 

/ A^ f (««)^) v 7 ?^ + / f r"^^ (<p £ ) e £ .(v £ ) v^dx £ 



JQ. C 



dx £ Vv £ ey £ 



(2.31) 



and 



/ g i i> e df<p e d E \fJgt<bf=[ P mii ' £ d £ y e %.{u £ )Jg^Ax £ V^G^. 

(2.32) 

From relation (2.28), it follows that the bilinear form associated with the left-hand side 
of ( l2~32l is ^-elliptic. 

Also for each h £ 6 V £ , the mapping 



P mij - £ d m \j/e £ llj (h £ ) v ^dx £ 



defines a linear functional on V P £ . Hence for each h £ £ V £ , there exists a unique T £ (h £ ) £ 
*F e such that 

/ <? ij : £ d £ T £ (h £ )d £ y £ Jg~tdx £ =f P mij ' £ d^e%, {h £ )Jg^dx £ Vy, £ eV £ 

Joe 1 JQf " J 
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(2.33) 



and that T £ : V £ — > I* 6 is continuous. 

In particular, it follows from d2.32t and the above equation that <p £ = T £ (u £ ) and 
eqs J2.3H and 12.321 become 



A^'^ ||/ (« £ )4 ; (v E ) v /^ck £ + / ^' e 4(T e (« e ))^ |y .(v e )^ck £ 

n e II NJ J ne IIJ 



f £ v £ v / g £ dx £ Vv £ eV £ , (2.34) 



/ g ,j ' £ d £ (T £ (u £ ))d £ y £ ^ £ ~dx £ =[ P mij ' £ d £ lV £ e £ ll .(u £ )J^dx £ 

Vy/ £ e^ £ . 



(2.35) 



Lemma 2.1. For each h £ G (L (O e )) , ?/;ere e^/ifi a unique G £ (h £ ) € V £ such that 



( A i ^ e e e k]]l (G e (h e ))^ lJ (v e )^dx e + [ ^ e d^T%G e (h e )))e^.(v e )^dx e 



h £ v £ ^/s £ <ix £ Vv £ eV £ 



(2.36) 



and that G £ : {L 2 (D. £ )y -> V £ ;i continuous. 

Proof. Let B £ (u £ ,v £ ) denotes the bilinear form associated with the left-hand side of 
eq. J2.34l >. Using J2.35i . we have 



B £ (u £ ,v £ ) = J^^ £ e £ k]]l (u £ )e £ y (v £ )^dx £ 

P"^- £ ^(r £ (« £ ))^ |/ (v E ) v ^ck £ 



A^' e 4 ||z ( M e )^.(v e )VF^ £ 



^^f(r £ (« £ ))5f(r £ (v E ))vFdx £ 



|n,(v e ) e ^.( M e )V^ck e 



+ / ^ £ «9, £ (r £ (y £ ))5f(r £ ( M £ ))^ £ dx £ 
B £ {v £ .u £ ). 



(2.37) 
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Also, using J2.35i and the relations 12.271 and J2.28l >. we have 



B £ {u £ ,u £ ) = / A' ikLr <- £ ..,",■>- 



^( M £ )4.(« e )^dx £ 

P'^d^iu^e^)^^ 



l!7W ' e "||| Z (" e )^|,-(" e )VFck e 



>C||w 



£||2 



lye- 



(2.38) 



Hence B £ (• ■ ■ ) is symmetric and V £ -elliptic. Hence by Lax-Milgram theorem, there exists 
a unique G e (h £ ) satisfying J2.36I . Letting v £ = G £ (h £ ) in J2.36t . we get 



v^4 ||z (G e (ft e ))4.(G e (/ i e )) v /^dx' 

F»y' e ^(r e (G E (A e )))4 ; .(G e (ft e ))v^ck £ 



/z £ G £ (/z £ )Vs £ d* £ 



Using J2.35i . it becomes 



^et^ih^e^m)^^ 



+ / £'J £ d £ {T £ {G £ (h £ )))d £ {T £ {G £ (h e ))) vV^ £ 



/z £ G £ (/! £ )v/pdx £ 



Using the relations J2.27I and J2.28i . we have 



||G £ (/ l £ )||2 £ <C £ ||G £ (/ ; £ )|| l/E ||^ll(^))3- 



Hence 



(2.39) 



(2.40) 



(2.41) 



\\G £ (h £ )\\ V e<C £ \\h £ \\ {L2m)i 

which implies that G £ is continuous. 



(2.42) 



It follows from J2.34i and the above lemma that u £ — G £ (f £ ). Since the inclusion 
(H^O 8 )) 3 ^ (L 2 (£2 £ )) 3 is compact, it follows that G £ : (L 2 (H £ )) 3 -> (L 2 (H £ )) 3 is com- 
pact. Also since the bilinear form B £ (- ■ ■ ) is symmetric, it follows that G £ is self-adjoint. 
Hence from the spectral theory of compact, self-adjoint operators, it follows that there 
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exists a sequence of eigenpairs (« m,e , <^ m,e )~ =1 such that 

/ n / ,,H ' e 4n/(« ffl ' e )ii. / (^)VFdx e 
+ ^p™'' £ ^(r £ ( M ffl ' e ))^ u .(v £ )v^ck £ 

= % m ' £ [ u m£ v £ J^dx £ Vv £ GV £ , (2.43) 

= j^P mii ' £ d £ n \if £ e £ y {u m ' £ )^Ax £ Vl^e^ e , (2.44) 

0<^ Lf <^- £ <.-.<^ £ < >oo, (2.45) 

j^ e u^ £ ^dx £ ^e 3 8 mn . (2.46) 

The sequence {m™' 6 } forms a complete orthonormal basis for (L 2 (Q.)) 3 . 
Define the Rayleigh quotient R(e)(v £ ) for v £ G V £ by 

= yA^ £ e m {v £ )e iy {v £ )^6x £ + yP^ 
(V) J aE vfvfVFdx £ 

Then 

|"«> e = m i n max fi £ (v £ ), 

W*€W«ve€W e \{0} 

where W £ denotes the collection of all m-dimensional subspaces of V £ 
3. The scaled problem 

We now perform a change of variable so that the domain no longer depends on e. With 

x = (x\,X2,xj) G £2, we associate x £ = (xi,X2,£xi>) G £2 £ . Let 

r = T6x(-i,i), r 1 = 7ix(-i,i), r ± = cox{±i}, 

r e = y e x(-l,l), r s = Ys>< (-1,1), 

r A , = r 1 ur+ur-, r eD = r+ur-ur e . 

With the functions Y p ' £ ,g £ ,A ijkl < £ ,P ijk < £ ^'J' 8 : £2 £ -» R, we associate the functions 
F> (e) , g £ , A''' w (e) , P^'* (e) , (e) : £2 -» R defined by 

P( £ )(x) :=P' £ (x £ ), g(e)(x) = /(* £ ), A"*(e)(jc) - A^ £ (x% 

(3.1) 

(e) (x) = P y *' e (x £ ) , ^ (e) (x) = £ ij ' £ (x £ ) . (3.2) 



(2.47) 
(2.48) 
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Assumption. We assume that the shell is a shallow shell, i.e. there exists a function 9 G 
C 3 (co) such that 

<j) £ (x h x 2 ) — (xi,x 2 ,£0(xi,x 2 )) forall (xi,x 2 ) G <a, (3.3) 

i.e., the curvature of the shell is of the order of the thickness of the shell. 
We make the following scalings on the eigensolutions. 

uT{x £ ) = e 2 u%(e) (*), v a (x £ ) = e 2 v a (x), (3.4) 

«™' £ (x £ ) - e«?(e)to, v 3 (x £ ) = ev 3M' ( 3 - 5 ) 

r £ (« m ' £ to)) - e 3 r(e)( M m (e)to), r £ (v(x £ )) = £ 3 r( £ )(vto), (3.6) 

§ OT ' e = e 2 ^ m (e). (3.7) 
With the tensors e% ■, we associate the tensors e.-ii . (e) through the relation 

1 1 1 y 1 1 J '"" 

4j(v e )(* e )=e 2 «i||/(e;v)(x). (3.8) 

We define the spaces 

y(i2) = {vG(// 1 (n)) 3 ,v|r =0}, (3.9) 

'P(£l) = {V/e// 1 (a),v/|r fD =0}. (3.10) 

We denote <p m (g) = r(e)(w m (e)). Then the variational equations (eqs ( !2.43> -( l246l ) 
become 

£1 

" P 3kI d, ( p m (e)e klll (e 1 v)^{ejdx 

a " 



+ e 



/ P ak '(e)d a( p m (e)e klll (e,v)^{ejdx 
Ja 



= $ m (e) / [e 2 <(e)v a + <(e)v3]v^(e)dx for all v e V(n). (3.11) 
^(e)d 3( p m (e)d 3¥ yMi)dx 



+ e / K 3a (e)(a a (p m (e)a3^ + ^<P m (e)^^)]vW)dx 
+ e 2 / ^(e)5 a (p m (e)^rV^)dx 



P M, ( e )5 3 ^||/(e," m (e))v / ^)dx 



e / [P ak '{e)d a \ l fe m {e,u m {e))]^{e)Ax for all ^€f(Q), (3.12) 

£2 



[£ 2 <(e)"«(e) + «3 (£)«§(£)] >/g(£)djc= 5,„„. (3.13) 
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4. Technical preliminaries 

The following two lemmas are crucial; they play an important role in the proof of the 
convergence of the scaled unknowns as e — > 0. In the sequel, we denote by C\ , C2 , • ■ ■ , C„ 
various constants whose values do not depend on £ but may depend on 0. 

Lemma 4.1. The functions e^j(e,v) defined in ( 13 ,8i are of the form 

e«||j3 (e;v) =e a p(v) + eV ce||j3 (e;v), (4.1) 

e a p(e;v) = l -{e ai {v)+E 2 el p {e-v)}, (4.2) 

e 3 ||3(e;v) = ^2e 3 3(v) ) (4.3) 

where 

e a p(v) = ^{d a vp + dpv a ) - v 3 d a p0 , (4.4) 

e a 3(v) = ^(<? a v 3 + <? 3 v a ), (4.5) 

e 33 (v) = ^v 3 (4.6) 
and there exists constant C\ such that 

sup max||«L(e;v)|| 0l n<Ci||v||i 1 o forallveV. (4.7) 

0<£<£o <*J 

Also there exist constants C2, C 3 and C4 such that 

sup max\g(x)- 1| <C 2 e 2 , (4.8) 
0<e<Eo xen 

sup max \A ijH (e)-A ijk '\ <C 3 e 2 , (4.9) 
0<e<eo 



where 



and 



A ijkl = X8 U 8 kl +il ^ik 8 jl + gUgjk) (4 10) 



A^MuMij > QMijMij (4.11) 



for < £ < £q and for all symmetric tensors (M, ; ). 
Proof. The proof is based on Lemma 4.1 of 1 2 1 . 



From relation i2.6\ and definition d3.2i . it follows that there exists a constant C5 such 
that for any vector (f,) £ M?, 

3 

#' 7 (e)^/>C 5 £r?. (4.12) 
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We assume that there exists functions P kl i and § l i such that 

sup max \P kij {e)-P kij \ <C 6 e, (4.13) 

0<£<£ X< = n 

sup max\£ ,ij (e)-£' ij \<C 1 s. (4.14) 



0<£<£o 



xeQ. 



Lemma 4.2. Let 6 G C 3 (o) be a given function and let the functions eij be defined as in 
J4.4t -( l4~6l . Then there exists a constant C% such that the following generalised Korn 's 
inequality holds: 

||v||i 1 n<C g |£||ey(v)||§ in | (4.15) 

for all v £ V(£l) where V(£l) is the space defined in ( I3.9I I. 
Proof. The proof is based on Lemma 4.2 of 1 2 1. 



5. A priori estimates 

In this section, we show that for each positive integer to, the scaled eigenvalues {<^ m (e)} 
are bounded uniformly with respect to £. 
Let <p € Hq(co). Then 



and 



Hence 



vcp := (-xidifp^-xidiq),^) e V(£2) (5.1) 
eap( v (p) = - x 3 d ap<P-<P d ap6, e a (v 9 )=0. (5.2) 



e a \\p(e,v<p) = -~x 3 d a p(p-(pd afi e + 0(e 2 ), (5.3) 

e a \\3(e,v<p) = 0(e), (5.4) 

e 3 p(e,v q> ) = 0. (5.5) 
We need the following lemma to prove the boundedness of the scaled eigenvalues. 
Lemma 5.1. There exists a constant C9 > such that 

|*(r(e)(vf))| 0| Q<q,|p|2,a., (5.6) 

\ed a (T(e)(v v ))\ ,n < C 9 \(p\ 2 ,m. (5.7) 



Proof. With the scalings J3.3i — (I3.7i . the variational equation (eq. (I2.33» posed on the 
domain Q. reads as follows: 
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For each h S (H l (Q.)) 3 , there exists a unique solution T{e){h) 6 (//'(H)) 3 such that 



/ <f 33 (e)a 3 7 , (e)(/!)^rVg(e)^ 
Jo 



+ e / [g a3 (s)(d a T(e)(h)d 3 yi/ + d 3 T(£)(h)d aW )}y^gJi)dx 

+ e 2 g«t (e)d a T(e) (h)d p wVgJijdx 

= [ P 3kl (e)d 3 Ye k]]l (e,h)^/gJ^dx 
Ja 

+ e [ P ak \£)daVe m {£,h)^&j&x Vy/e*F. (5.8) 
Jo. 

Taking h — and y/ = T(e)(v,p) in the above equation, we have 

J^ 3 (e)d 3 T(e)(v (p )d 3 T(e)(v (p )^gJi)dx 

+ eJ^ a \e){d a T{e){v <s> )d 3 T{e){v <s> ) 

+ d 3 T(e)(v <l> )d a T(e)(v (p ))]^gJi)dx 

+ e 2 ( ^(e)d a T(e)( V(p )dpT(e)( V(p )^Mejdx 



= j^ kl (e)d 3 T{e){v (p )e kV (e 1 v ip )^(e)Ax 



+ e / P akl (e)d a T(e)( V(p )e k]]l (e,v <p )^)dx. (5.9) 



Using the relations J4.12> and J5.2I >— d53i . it follows that there exists a constant Cg > 
such that 

\MT(£)(vM,n + \£MT(£)(vM.n 

< Cj{|a J r(e)(v,)| ,n|p|2,o» + \ed a T(e)(y v )\op\<p\ 2t0 } (5.10) 

and hence the result follows. 

Theorem 5.2. For each positive integer m, there exists a constant C{m) > such that 

r(e)<C(m). (5.11) 

Proof. Since problem d3 .111 was derived from ( 12.431 after a change of scale, we still 
have the variational characterization of the scaled eigenvalues t, m {e). Let V m denote the 
collection of all ra-dimensional subspaces of V(£l). Then 

^ m (e) = min max ^7^4 > (5-12) 
b v ' wev„, vew D(e)(v,v) 
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where 

N(e)(v,v)= [ A^%| / (e,v) e!l | 7 (£,v) v ^(i)dx 

+ / P 3kl d 3 T(e)(v)e klll (e,v)yMe)dx 
Ja 



+ e / P akl d a T(e)(v)e klll (e,v)^{e)dx, (5.13) 
a 



D(e)(v,v)= / [e 2 y a v a + v 3 v3]V^(e)ck. (5.14) 



a 

Let W m be the collection of all m-dimensional subspaces of Hq((o). Let W € W m . Define 
W={v 9 |(pGW}. (5.15) 
It follows that W e V m . Hence, it follows from (15.121 that 

< min max ^ E j ( / y ' Vy j . (5.16) 
wew m <peW'D(£)(v (f ,,v ¥ ,) 

Now, 



D(e)(v 9 ,v 9 )= / £j [e 2 ^|5 a <p| 2 + |<P| 2 ]V^(e)^ 



> / <p 2 dco. (5.17) 



(5.18) 



Using the relations J5.3i — < l5T5t and Lemma 5.1, it follows that 

A ,; '%| / (e,v v )e !l | i (e,v ( p)V5(e)ck<C / |Ap| 2 da), 

i J (a 



^P 3k %T(e)(v v )e m {e,v< p )^f^dx < C / |A<p| 2 dco, (5.19) 



e / P aW a„r(e)(v (| ,)^|| / (e,v ( p) v ^(e)dx<C / |A<p| 2 dco. (5.20) 

Hence 

§«(«)< Cidn maA |A f d<a 

< CA m , (5.21) 
where A" ! is the mth eigenvalue of the two-dimensional elliptic eigenvalue problem 
A 2 u = Am in co 

U = d v u = on dco. (5.22) 
This completes the proof of the theorem on setting C(m) = CX m . 



346 NSabu 
6. The limit problem 

Theorem 6.1 . (a) For each positive integer m, there exists u m € H l (Q.) , <p m e L 2 (£1) and 
E, m et such that 

u m (e)^u m inH l {Q.), <p m (e)^<p m inL 2 (£2), (6.1) 

(ed l( p' n (e),ed 2( p m (e),d 3( p m (e)) - (0,0,^<p m ) inL 2 (£2), (6.2) 

§l»( £ ) (63) 

(b) Define the spaces 

V H (co) = {(ri a )e(H l (co)) 2 ;ri a =OonY Q }, (6.4) 

V 3 (fi)) = {773 e H 2 (co);ri 3 = d v i} 3 = on y }, (6.5) 

VjfL = {v e H\£l)\v = rj a -x 3 ^773, (77,) e V H (co) x v 3 (©)}, (6.6) 

^ l = {weL 2 (a),d 3 xi/eL 2 (Q.)}, (6.7) 

<F /0 = {y € L 2 (Q.),di W e L 2 (£2), v/|r± = 0}. (6.8) 

Then there exists £™) € V// x V 3 (co) smc/z that 

u^ = Q-x 3 d a ^ and «3 n = C 3 m , (6.9) 
3a/3 

9" = (1-^)^-5^3" (6.10) 
and (£ m ,§ m ) € Vg x V3 x R wiris/ies 

= <? m / C 3 m T?3dfi) V773 Gy 3 (co), (6.11) 

J^n e ap d p r la d(O = Vri a eV H (co), (6.12) 



vv/zere 



"«/»(£) = - { 3^+^) A ^ 5 -P + T a »^} (6 ' 13) 

n%(Q = ^l^(OS a p +4^/3(0 (6-14) 

33 = I j p3a3p3a3 + _!_ p 333 p 313 ^33 (fi _ 15) 

A +2ju 

3«/3 = p 3«/3 _ ^^p333 5 a/3 (fi lfi) 
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Proof. For the sake of clarity, the proof is divided into several steps. 
Step (i). Define the vector and the tensor K m (e) = {K™(e)) by 

$*(£) - {£d l(p m {e),ed 2(p m {e)A(p m {e)), (6.17) 

K^(e)=e ap (u m (e)), KUe) = \e ai {u m (£)), K? 3 (e) = ±e 33 (u m (e)). 

(6.18) 

Then there exists a constant Cio > such that 

ll" m (e)lli.n<C 10 , |^(e)|o,n<C 10 , |$"(e)|o,n < Cio (6.19) 

for all < e < £o- 

Letting v = m™(£) in ( 13.1 11 1, we have 

/ A^( e K|| / ( e )( M '"( £ )) ei1 | 7 ( e )( M '"( e )) v ^(i)dx 

Ja " " 

+ / P 3H ( £ )^<P m (eK||/(e)("'"(e))v^(i)dx 

+ e / P akl (e)d a < P m (e)e m (e)(u m (e))^(i)dx 
Ja 

= <?'"(£)/ [e 2 «™(e)«™(e) + <(£)<(£)] V^)ck. (6.20) 

Ja 

Letting i/a = (p m (e) in ( 13 . 1 21 and using it in the above equation, we get 

/ A-j k '(e)e klll (ey n (e))e illj (e, u m (e))^{i)d x 
Ja 

+ J^{e)^{e)ipJ{e)^ME)Ax 
= r(e)/ [e 2 <( e ).<( e ) + <(eK(e)]^(i7dx. (6.21) 

Using the coerciveness properties ( 14. 1 1 i and J4.12> . the inequality (a — b) 2 > a 2 /2 — b 2 
and the generalized Korn's inequality (14. 151 . we have for e < minjeo, 1}, 

/ A'i k \e)e kV {e,u m {e))e iy {e.u m {e))^{e)d x 
Ja 

+ jf ^(e)C(e)9f (e)V&)dx 
>C 11 Elk I ||>(e I » m (e))llo,n + Ci 1 EllC(e)llo,n 
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--CnZ\\e a p(u"\e)) + e 2 el p (e,u m (e))\\l a 
a,j8 



+ 2C„£ 

+ C n 



l -e a2 {ir{e)) + ee^{e,u m {e)) 



2 

o.n 



e 33 (« m (e)) 



-CiiEHC(e)llo,a 



> c n|^Ll^"( e )laa-Cf(2 £ 2 + e 4 )|| M m ( e )||L| 
+CiiEl|^T(e)||§ |0 

i 

^ C 'i|^Lll^(» ffl (e))llo,n-3 e 2 C?||« ffl ( £ )||2 a | 

+ C 1 iEH^"(e)||^ 
>C 11 {i(C 8 )- 2 -3 e 2 C 2 |||« m (e)llm + CiiEllC(e)llL- 



(6.22) 



Combining eqs ( 16. 2H and J6.22I with relations < I3 . 1 31 and ( 15.1 It . we get the relation 
( EH9l l. 

5fe/7 (ii). From Step (i) it follows that there exists a subsequence (<p" ! (£)) and (<p™) 6 
L 2 (f2) such that 

(e^ 1 (p" 1 (e),ea 2 (p" 1 (e),53(p m (e)) - in (i 2 (£2)) 3 . (6.23) 

Since r e £) contains r~, we have 



(p m (e){x l ,x 2 .x 3 ) = / a 3 ^ m (e)(x 1 ,x 2 ,i)ds 



(6.24) 



and it follows that ||(p m (e)||o,n < V2||^3<P m (e)||o.a- This implies that <p m (£) is bounded 
in L 2 (Q.). Therefore there exists a <p m in L 2 (£2) and a subsequence, still indexed by e, 
such that <p" ? (e) converges weakly to (p m . Hence it follows from J6.23b that 



(ed l( p m (e),ed2<p m (e),d 3 <p m (e)) - (0,0,5 3 <p m ). 



(6.25) 



Sfe/} (iii). From Step (i) it follows that there exists a subsequence, indexed by e for 
notational convenience, and functions u m £ V (£2) and A™ e (i 2 (£2)) 9 such that 



u m (e)^u m in H\D.), K m {e)^k m inL 2 (£2)ase^0. 



(6.26) 



Then there exist functions e ff 1 (fl)) and £™ G H 2 (») satisfying £/" = <9 V C 3 " = 
on 7o such that 



u% = Q-x 3 d a t;? and < = 



(6.27) 
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and 

V m z l,, m '\ V m r>3oi3 ~\ m 

K ap= e ap{u ), K a3 = -—P d 3 (p , 

*™ = -TT2t {pmd3(pm + XR PP ] - (6 ' 28) 

From definition ( I6.18> and the boundedness of (Rfj(e)), we deduce that 

lka3(M m (e))llo,n<eCi 3 and ||e 3 3(" m (e))||o,n < e 2 Ci 3 , 
where e i; (v) = h(diVj + djVj). Since norm is a weakly lower semicontinuous function 

IM«™)lta < liminf \\e i3 {u m (e)\\ ,n = 0, (6.29) 

we obtain e,3 (u m ) = 0. Then it is a standard argument that the components u'" of the limit 
are of the form (I6.27i . 

Since u m {e) — 1 u m in H l (Q.), definition ( I4.4t of the functions e a p(v) shows that the 
function (e) = e a j3(w m (e)) converges weakly inL 2 (£2) to the function e a p(u m ). 
We next note the following result. Let w £ L 2 (£2) be given; then 

wd 3 vdx = for all v E H l (£2) with v = on T , then w = 0. (6.30) 

! 

Multiplying (13 . 1 1 i by £ 2 , taking (v a ) = and letting e — > 0, we get 

(X%* a + (X+2^)K 33 +P 333 d 3 (p m )d 3 v 3 dx = (6.31) 

! 

which implies {XK% C + (A + 2fi)K 33 + P m d 3 (p'") = and hence the third relation in 

d08l follows. 

Again, multiplying J3 .111 by e, taking V3 = and letting e — > 0, we get 

(M^ 3 + P 3a3 d 3 <P m )^v'ad* = (6.32) 

which implies {llK% 3 + P 3a3 d 3 (p m ) = and hence the second relation in (I6.28> follows. 

Step (iv). The function <p™ is of the form J6. 10i . 
Letting e — > in eq. (13. 121 . we get 

{P 3a f i K™ l3 -<g' 33 d 3 <p m )d 3 \i/dx = Vy/E^iO.). (6.33) 

Since D(£2) is dense in 'P/o (and hence in *P(£2)) for the norm ||.||<p,> eq. (I6.33> is 
equivalent to 

d 3 (P 3ap K" p - S 33 d 3 (p m ) = in D'(O) (6.34) 
which implies that {P 3a ^K^ p ~ <? 33 <? 3 <p m ) = d\ with c/ 1 e D(o)). Then 

^<P'" = ^-MD -^d aP Q"} - (6.35) 
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which gives 

«p»> = P^- [x3 e ap{n - x \d a ^\ - ^Ld l +d°. (6.36) 
Since (p m satisfies the boundary conditions q>^ + = <pj^_ = 0, we have 

d° = ^d ap ^, d l =p^~e afi {n- (637) 
Thus the conclusion follows. 

Step (v). The function (£•") satisfies (I6.1H and J6. 12i . 
Taking v € Vkl an d letting e — > in ( 13. 1U we get 

A»^£ aj8 (v)d* + / P ia Pd,(p'"K ap (v)dx = / u n i ■ v 3 dx. (6.38) 



n Jo. Ja 

Replacing u m and K'j" by the expressions obtained in (I6.27> and (I6.28> . and taking v of the 
form 

v« = r/a -X3d a r]3 and v 3 = tj 3 
with (tj,-) € V//(co) x Vi((o), it is verified that (16.381 coincides with eqs ( I6.11l l and ( I6.12> . 

Sfep (vi). The convergences u m (s) — 1 m" 1 in i/ 1 (£1) and <p™(£) — »■ <p m in L 2 (£l) are strong. 

To show that the family (w m (e)) converges strongly to u m in H l (fi), by Lemma 4.2, it 
is enough to show that 

£,■>"'(£)) ^e, 7 (" m ) inL 2 (£2). (6.39) 

Since e^u" 1 ) = and 

^\\e l] {u"\e))-~e lj (u m )\\l ) n 

= L ll^(e)-^llo,a + 2e 2 Lll^ye)||ia + e 4 |l^33(e)llaa, (6-40) 

a,p a 

convergence ( 16.391 1 is equivalent to showing that 

K' n (e)^K m inL 2 (£l). (6.41) 

We define a norm on (L 2 (£2)) 9 x (L 2 (£2)) 3 by letting for any matrixM £ (L 2 (£l)) 9 and 
any vector % € (L 2 (fi)) 3 , 



, , ( ^ 1/2 

A'^'m : Mv^(i)dx + J g' J XiXjVg(£jdx \ . (6.42) 



LetX m (e) be the norm of (K m (e), edi (p m (e),e d 2 (p m (e),d 3 (p'"(e)) in ( L 2 (H) ) 12 . Using 
the weak convergence equation (eqs (I6.25> and ( I6.26t ) and the relation ( I6.28t . it can be 
shown that 

\imX'"(e)=X m = ( [ A ijkl K m :K m dx+ [ g 33 {dj,<p m fdx \ ' (6.43) 
£^0 \Jn Jn I 
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which is the norm of (AT"', 0,0, d^(p m ). Since we have already proved that (K m (e), 
e<9i<p'"(e), ed 2 (p m {e), d 3 (p m (e)) converges weakly to (£,0,0,<? 3 <p m ) in (L 2 (H)) 12 , we 
have the following strong convergences: 



K m (e) -> K m strongly in (L 2 (£2)) 9 , 

( £ d l <p m ( £ ), £ d 2 <p"'( £ ),d 3 q> m (e)) (0,0,5 3 (p m ) strongly in (L 2 (£2)) 3 



(6.44) 



(6.45) 



Hence u m { £ ) converges strongly to u m in H l (Cl) and since (p m ( £ ) — (p m is in ^P/o, 
the equivalence of norms || y||<p, and — > | c*3 1//" | n in ^P/o proves that (p'"( £ ) converges 
strongly to (p m in L 2 (H). 



Equation ( I6.12t can be written as 



dpT} a dC0 



2X/x 
_X+2n 

Clearly, the bilinear form 



dpTjadCO. (6.46) 



2AjU 
A +2ji 

2Xjx 



e PP (QS a p + 2/j,e a p(Q 



dp7] a d(0 



is Vh(co) elliptic. Also for a given £3 G V 3 (co), the functional 



dco (6.47) 



j^-{da6d„k)8 a p+li(d a ediik + di l 6d a &) 



dpriadeo 
(6.48) 



is continous on Vh{co). Thus, given £3 € V^ffl), there exists a unique vector (£ a ) € V/y (ft)) 
such that 



(6.49) 



We denote by e V g (o) x V 3 (ft)) the vector (£«, £5). In particular, r£™ = (C« , C 3 m )- 
Substituting this in (I6.1H . we get 



b(Q\ 173) = 4 m / C m T} 3 d(B forallT}3Gy 3 (a)), 



(6.50) 



where 



*(£»,1fc) 



2 f p^PplP* 



3 id) p 



33 



(6.51) 
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Lemma 6.2. The bilinear form b(- ■ ■) defined by (16.5 It is VH(co)-elliptic and symmetric. 
Proof. It follows from Lemma 6.2 in 1 8 1 that the bilinear form &(•••) defined by 

&(6,T?3) = - / m ap (t; 3 )d a pT l3 dcQ+ ( n e a JT^ 3 )d a p9T] 3 da) (6.52) 

is V//(o)) -elliptic and symmetric. Hence it is clear that &(•••) is also Vh{(£>) -elliptic and 
symmetric. 

Lemma 6.3. Lef (C,™ ,E, m ),m > 1, fee f/ie eigensolutions of problem ( 16.511 1 found as limits 
of the subsequence (u m (e) (e)) ,m > 1 of eigensolutions of the problem (13. lit . 77zen 
f/ze sequence (§ m )~ =1 comprises all the eigenvalues, counting multiplicities, of problem 
( 16. 5U and the associated sequence (C™)m=i of eigenf unctions forms a complete orthonor- 
mal set in the space V 3 {(o). 

Proof. The proof is similar to the proof of Lemma 5.4 in 0. 
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